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Probability distribution function for a solid with vacancies 
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Expression for probability distribution is got taking into account a presence and removal of de- 
generacy on the microstates. Its application allows to describe the process of melting of solids, as 
saltatory phase transition of the first kind without bringing in of concept of the order parameter. 
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For the case of a solid with vacancies the expressions 
for the thermodynamic probability W and the configura- 
tion entropy S c in supposition of the total degeneracy of 
the microstates had formerly got [l[ 



W = 



(N + n)\ 
NW. ' 



S c = k B lnW, 



(1) 



(2) 



where N is a total number of atoms in a crystal, n is a 
number of vacancies, k B is the Boltzmann constant. 

The total probability of the state with n vacancies 
contains a limiting exponential multiplier, probability of 
Gibbs for this microscopic state [H- 



/w=!« P (-fM), 

Z ks-l 



(3) 



where Z is the normalized statistical sum, U(n) is the in- 
ternal energy taking into account existence of vacancies, 
T is the temperature. In the case of total degeneracy 
the internal energy is the linear function of number of 
vacancies [l[ 



U (n) = Uq + u n, 



(4) 



where Uq is the internal energy without taking in account 
of vacancy contribution, uq is the energy of a vacancy 
which in this case is a constant for all possible configu- 
rations of atoms of a system. 

The condition of total degeneracy is executed more or 
less exactly for the low concentration (number) of vacan- 
cies, when it is possible to ignore the interaction between 
vacancies. At the same time, this condition is not ex- 
ecuted for those configurations for which vacancies are 
located close to each other or merges at all (bi- vacancies, 
triple vacancies, vacancy pores) and at the high concen- 
tration or density of vacancies. In these cases it is neces- 
sary to take into account the removal of condition of total 
degeneracy due to the interaction between vacancies. As 
it is known, coupling interaction reduces the total energy 



of a system that reduces the effect of action of limiting 
exponential multiplier in Eq. [3] and resulting to possi- 
bility of appearance of long-living configurations, which, 
thus, can resist more numerous configurations with more 
high energy. 

The probability distribution function (PDF) with total 
(internal) energy E\ equals [|| 



fin) 



w{Ei) 



exp(- 



Ei 
k B T 



), 



(5) 



where w(E{) is distribution of states with energy or num- 
ber of the microstates (configurations) with energy of Ei, 
Z is statsum over all energetic states of the system 



Allstate 



Z= J2 w(Ei)exp(- 



i=i 



k B T> 



(6) 



The states are numbered in order of its energy growth 
of Ei+i > Ei. Distribution w(E{) depends of the number 
of vacancies and of symmetry of their location (ordering). 
On this account, it is impossible to fix one-to-one corre- 
spondence between internal energy and number or con- 
centration of vacancies. Finding of function w(E{) is an 
intricate combination problem, determination of power 
spectrum of Ei, as the set of acceptable energies for I 
state, is also a stubborn problem. However, such one-to- 
one correspondence can be got for the average value of 
internal energy U on all of the states at the fixed num- 
ber of vacancies n. Using expression for PDF <(5j) , it is 
possible to write down (see closed idea in Q, p. 147) 



U(n) = 



All states (n) 

- y 

z ^ 

1=1 



E lW {Ei)e^{--^=), (7) 
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Here unlike a case of Eqs. [T]-|4]the internal energy is not a 
linear function, but is the function of general view. The 
increase of fraction of the symmetric low-energy states 
diminishes the average value of the internal energy. The 
transition to the states, at which vacancies at the high de- 
gree of their concentration can pass to higher-energy ex- 
tended (and more mobile) states [5|-l7|, results in growth 
of the internal energy. For the reflection of all of these 
properties of the internal energy it presents in a form of 



2 



polynomial with alternating signs 

K r-nfc-i 
U = U + J2—^ — u k-m\ (8) 

k=l 

Ignoring degeneracy within of one set of the states with 
the same number of vacancies n, for such average internal 
energy we can write down the effective function of distri- 
bution of probabilities in a form Eq. [3] with determina- 
tions ((TJ, @ and flSJ. Using such PDF one can describe 
processes of melting of solids, as phase transitions of the 
first kind 0. 

Values n, at which PDF has extreme values (equilib- 
rium states), can be found from the follow transcendent 
equation || 

fc=l k=l 

where 

u ^^r = E(-i) w ( 10 ) 

fc=0 

is energy or chemical potential of a vacancy. 

The first two terms in Eq. [§] are sums S of the a 
slowly divergent harmonic series. This part of the equa- 
tion depends only on system size and doesn't depend on 
material parameters. It is the fundamental part of Eq. 
[5] decreasing with the growing parameter n (curve 1 of 
Fig. rjJ[21E])- The last term in Eq. [^depends on material 
parameters through coefficients Uk- 

Now let us consider a problem in varying degree of 
approximation. In the linear approximation on the num- 
ber of vacancies for the internal energy or at constant 
energy of vacancy there is only one solution of Eq. R5J 
which coincides with the well known Frenkel solution [l| 
(see Fig. [T|). At lowering of the energy of vacancy or 
at the increase of temperature the solution is continu- 
ously displaced to region of higher number (or concen- 
trations) of vacancies. We can formally get any value of 
vacancy concentration, but, it is clearly, that too high 
vacancy concentrations about 100 percent lose physical 
sense and theory becomes useless. However Frenkel had 
proposed to suggest not too high vacancy concentrations 
about 10 percent, as transition to the liquid state. At the 
same time, the solution offered by him (in approximation 
of independence of energy of vacancy from the vacancy 
concentration) can not explain or describe the process of 
melting of solid, as saltatory phase transition of the first 
kind 

In the quadratic approximation on the vacancy num- 
ber for the internal energy or in the linear approximation 
for energy of vacancy Eq. [9]has already two solutions in a 
region of high energy of vacancy or in region of low tem- 
peratures (see Fig. [5]). One of them in region of lower 
concentration of vacancies, actually, coincides with the 
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FIG. 1. Solution of Eq. [5]at constant energy of vacancy. The 
fundamental curve S (1), the energy of vacancy u (2) and 
PDF (3): u = 0.1; 0.07; 0M(eV). Other parameters equals 



,T = 


300K. 


The right-hand scale is for PDF. 
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FIG. 2. Solution of Eq. [9] at linear dependence of va- 
cancy energy of the vacancy concentration. The fundamen- 
tal curve S (1), the energy of vacancy u (2) and PDF (3): 
u = 0.12; 0.105; 0.1; 0.085(eV), ul = 0.0002(eF). Other pa- 
rameters equals zero, T — 300K. The right-hand scale is for 
PDF. 



Frenkel solution, but it is some displaced due to inter- 
action between vacancies. The second solution in region 
of higher concentration of vacancies describes the equi- 
librium or stationary solution of problem also; however, 
it corresponds to a minimum of PDF and is unsteady. 
Probability of realization of this state is less not only in 
relation to the stable stationary state but also in rela- 
tion to any non-equilibrium state. In addition, for the 
low energy of vacancy or for high temperatures we get 
in a region in which the solutions of Eq. [§] are absent at 
all. The last circumstance testifies unsuitability of the 
approximation in this region, and it is need to resort to 
a higher approximation. 

In the cubic approximation on the number of vacancies 
for internal energy or in the quadratic approximation for 
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FIG. 3. Solution of Eq. [5] at quadratic dependence of energy 
of vacancy from the vacancy concentration. The fundamental 
curve 5* (1), energy of vacancy u (2) and PDF (3): uo : 
0.135; 0.13; 0.125; 0.115(eV). The rest parameters are m = 
0.6 ■ 10- 3 (eV), u 2 = 0.125 ■ nT 6 (eF), T = 300K. The right- 
hand scale is for PDF. 



to the nonlinear contributions. At lowering of the energy 
of vacancy or with growth of temperature equation will 
have three solutions, one of which (left) is the modified 
Frenkel solution, second one (intermediate) is unsteady, 
and third one (right) in region of high concentration of 
vacancies can be possible to examine as proper the liquid 
state of a matter. In this case we have an equilibrium 
coexistence of solid and liquid phases of the matter. At 
the yet greater increase of temperature we pass to the re- 
gion of existence of only one solution which corresponds 
to clear-melted matter. Thus a transition from the solid 
state to liquid one is carried out as a phase transition of 
the first kind. 

In the article the deduction of a probability distribu- 
tion function is considered for the case of general de- 
pendence of the internal energy from the number (con- 
centration) of vacancies. Within the framework of such 
generalization it is shown the melting of solids can be de- 
scribed, as a saltatory phase transition of the first kind 
without bringing in of a concept of order parameter. 



the energy of vacancy Eq. [§] can have one or three solu- 
tions already (see Fig [3]). In region of high values of the 
energy of vacancy or low temperatures equation there is 
one solution, which is the Frenkel solution modified due 
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